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ABSTRACT 

We present the results of SPITZER mid-infrared spectroscopic observations 
of two highly-obscured massive X-ray binaries: IGR J16318-4848 and GX 301- 
2. Our observations reveal for the first time the extremely rich mid-infrared 
environments of this type of source, including multiple continuum emission com- 
ponents (a hot component with T > 700 K and a warm component with T ~ 
180 K) with apparent silicate absorption features, numerous H I recombination 
lines, many forbidden ionic lines of low ionization potentials, and pure rotational 
H2 lines. This indicates that both sources have hot and warm circumstellar 
dust, ionized stellar winds, extended low-density ionized regions, and photo- 
dissociated regions. It appears difficult to attribute the total optical extinction 
of both sources to the hot and warm dust components, which suggests that there 
could be an otherwise observable colder dust component responsible for the most 
of the optical extinction and silicate absorption features. The observed mid- 
infrared spectra are similar to those from Luminous Blue Variables, indicating 
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that the highly-obscured massive X-ray binaries may represent a previously un- 
known evolutionary phase of X-ray binaries with early-type optical companions. 
Our results highlight the importance and utility of mid-infrared spectroscopy to 
investigate highly-obscured X-ray binaries. 

Subject headings: circumstellar matter — infrared: stars — stars: emission-line, 
Be — X-rays: binaries — X-rays: individual (GX 301-2, IGR J16318-4848) 



Introduction 



Recently, a large number of highly-obscured (e.g., A^h > 10^'^ cm ^) massi ve X-ray bina- 



ries h ave been discovered by the INTEGRAL hard X-ray (> 15 keV) satellite (jWinkler et al. 



20031 ). The prototypical case is IGR J16318-48 48 which shows variable high obscuration i n 
the X-ray, sometimes reaching A^h — 2 x 10^^ (ICourvoisier et al.ll2003l : IWalter et al.ll2003l ). 
Its bright optical and ne ar-infrared (IR) counterp art is an early B-type supergiant star with 
numerous emission lines (IFilliatre fc Chatyll2004j ). Interestingly the obscuration toward IGR 
J16318-4848 obtained in the optical and near-IR wavebands {Ay ~ 18) is almost two orders 
of magnitude smaller than that inferred from the X-rays, which suggests that the extreme 
obscuration seen in the X-ray is intrinsic only to the X-ray source. However, Ay ~ 18 is still 
greater than the interstellar obscuration, indicating the existence of substantial circumstellar 
material around the supergiant companion. 

In order to fully understand the implications of the INTEGRAL discoveries — specif- 
ically if they imply existence of a separate class of highly-obscured X-ray binaries — we 
must investigate any similarities b etween the new INTEGR AL sources and previously known 
sources. As already suggested by iRevnivtsev et al.l (120031 ) . the X-ray pulsar GX 301-2 ap- 
pears similar to the new INTEGRAL highly-obscured massive X-ray binaries: it has variable 
X-ray obscuration of A^h — 10^^-10^^ cm~^, the compact sour ce is a neutron sta r, and the 
optical companion is an early B-type supergiant ( o r hyp ergiant; iKaper et al.lll995l ). like IGR 



J16318-4848. Recently, iKaplan. Moon, fc ReachI (|2006| . ; hereafter Paper I), using optical 



near-, and mid-IR (< 20 /im) spectral energy distributions (SEDs), have found that both 
sources have strong mid-IR excesses that they identify as continuum emission from hot dust. 
This is suggestive that they have very similar circumstellar material which may be related 
to their strong X-ray obscuration. In this Letter, we present the results of SPITZER mid-IR 
spectroscopic observations of IGR J16318-4848 and GX 301-2, showing that both sources 
indeed have very similar rich mid-IR properties previously unknown for X-ray binaries. 
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2. Observations and Data Reduction 

We ob served IGR J16318-4848 and GX 301-2 with InfraRed Spectrogrpah (IRS) aboard 



SPITZER flHouck et al.l 12004 : IWerner et al.l |200J) on 2005 September 8 and 2005 July 2 



respectively. The duration of the total observations was 1.32 (IGR J16318-4848) and 0.82 
(GX 301-2) hours with the all six IRS modules in operation in the standard staring mode, 
obtaining two nods of spectra separated by 1/3 of the slit length. For the basic data reduc- 
tion, we used the co-added Basic Calibration Data (BCD) produced through the standard 
SPITZER Science Center data reduction pipeline (S13.2.0), and also the SPITZER IRS 
Custom Extraction (i.e., SPICE) software (vl.3-betal). In the low-resolution mode where 
the slit lengths are large (i.e., 57" for 5.2-14.5 fim and 168" for 14.0-38.0 /xm), we sub- 
tracted out background emission of one nod using the spectrum of the another nod, and 
combined the two background-subtracted spectra for the final spectrum. However, in the 
high-resolution mode, the small slit length (i.e., 11"3 for 9.9-19.6 fim and 22"3 for 18.7- 
37.2 fim), together with the high brightnesses of the sources, made it impossible to subtract 
out the background. Therefore, we just combined the spectra from the two nods for the 
final high-resolution spectrum without any background subtraction. In addition, we used 
IRSCLEAN (vl.5) to mask and clean rogue pixels from the co-added BCD images, and then 
interpolated over neighboring normal pixels. 



3. Low-Resolution Spectrum: Continuum and Broad Features 

Some emission features, such as ionic forbidden lines, H2 lines, and Polycyclic Aro- 
matic Hydrocarbon (PAH) emission, present in the low-resolution spectra appear extended. 
Here we consider just the emission from point sources, and defer discussion of the extended 
emission to a separate paper. Figure [T] presents the background-subtracted low-resolution 
spectrum of IGR J163 18-4848, as well as the best-fit spectra of the continuum emission and 
the residual of the SED fit (see below). Note that the low- resolution spectrum is extracted 
for the continuum source (i.e., IGR J16318-4848) based on the standard point source extrac- 
tion procedure of SPICE. Several features are prominent in Figure [T|, including complicated 
multi-component continuum emission, broad and strong absorption around 9.7 and 18 fim, 
and numerous emission lines. The spectrum of GX 301-2 (Figure [2]) is similar, with strong 
continuum emission, silicate absorptions feature around 9.7 /im, and several emission lines. 

As in Figures [1] and [21 the silicate absorption features and numerous emission lines 
make it difficult to carry out any meaningful independent SED fits of the shorter wavelength 
(i.e., A < 8 /im) data of the IRS spectra. Therefore, instead, we first examined if the IRS 
SEDs at the shorter wavelengths are consistent with what have been found in Paper I using 
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only the line-free regions. According to Paper I, the SEDs of IGR J16318-4848 and GX 
301-2 at A < 8 yum have contributions from hot (i.e.,T ^ 1040 K for IGR J16318-4848 and 
720 K for GX 301-2) dust emission and stellar emission of a BOI star with temperature of 
26000 K, while the contribution of free-free emission is negligible. Our IRS SEDs of the 
line-free regions in Figures [1] and [2] at A < 8/im were fitted very well with those parameters, 
resulting in the reduced chi-square xl < 1- We used the extinction Ay = 18.5 and 7.1 for IGR 
J16318-4848 and GX 301-2, respectively, as in Paper I. However, in the longer wavelength 
(i.e., A > 22 /im) range, we were unable to fit the IRS SEDs with the given hot dust and 
stellar parameters because the best fit resulted in unacceptably large chi-square values. We 
therefore performed SED fits with three components (i.e., hot dust, stellar, and warm dust 
component) using the entire line-free SEDs in the A < 8 /im and > 22 /im range. Over the 
fits the parameters of the warm dust component were treated as free parameters, while those 
of the hot dust and stellar component were fixed to be the values in Paper I. As a result, we 
obtained the warm dust temperature of 190 K (IGR J16318-4848) and 170 K (GX 301-2) 
with xl - 1-2 (IGR J16318-4848) and 1.5 (GX 301-2). We used the extinction coefficients 



compiled in iMathis I (120001 ) excluding the 8-22 /im range in order to avoid the uncertainties 
of the coefficients associated with the silicate absorption features. We instead interpolated 
the extinction coefficients at A < 8 and > 22 /im to calculate the silicate absorption-free 
extinction coefficients in the A = 8-22 /im range, and used them in the fits. The main 
panel of Figures [1] and [2] presents the results of our best SED fit obtained by excluding the 
silicate absorption features, where we can confirm that the fit matches nicely the observed 
SED except the A = 8-22 /im range. In the bottom panel, on the other hand, we show 
the redisuals of the fit obtained by subtracting the fiuxes computed by the best fit from the 
observed values. For this residual construction, in order to see if our best SED fit is consistent 
with the observed silicate absorption features, we combined the best-fit parameters obtained 



at A < 8 and > 22 /im with the extinction coefficients in the A = 8-22 /im range (IMathis 



20001 ). and calculated the silicate absorption feature-associated fiuxes expected by the best 
SED fit in the A = 8-22 /im range. The residuals in Figures [1] and [2] show that overall 
our best SED fit agrees with the silicate absorption features reasonably well, although there 
appears to be some small discrepancies for IGR J16318-4848. 

The masses of the d ust components were estimated under the optically thin assumption 



(IHildebrand et al.l 119771 ): Mdust = Fxd / KxBx{Tdust)i where d is the distance to the source. 



K\ is the dust mass absorption coefficient, and Bx(T) is the Planck fun ction at temperature 



T. Using the interstellar mass absorption coefficients of iDraind (120031 ) at 8 and 25 /im, the 



hot and warm dust masses are estimated to be M^ot — 4.5 x 10 ^ c?| Mq and M^a^^ — 3.1 
X 10"^ dl Mq for IGR J16318-4848; M^ot ^ 5.7 x 10"^ dj Mq and M^arm ^ 1-4 x 10"^ dj 
Mq for GX 301-2, where d^ is the distance to the sources normalized by 5 kpc. 
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The silicate absorption features seen in Figures [T] and [2] have a tight hnear correlation 
with the optical extinction in the diffuse interstellar medium (ISM) , although the correlation 
has been reported to break down in dense enviornment where the si licate feature doe s not 
show a monotonic increase with extinction at Ay > 12 mag (e.g., IChiar et al.l 120071. and 
references therein). The correlation in the diffuse ISM is Ay/rgj = 18.5 (jDraind l2003l . and 
references therein; rg.y is the optical depth of the silicate absorption at 9.7 /im). If we adopt 
the above correlation given in the diffuse ISM, the intensity ratio of the observed spectrum to 
the unabsorbed spectrum in Figures [T] and [2] at 9.7 /im correspond to Ay ~ 18.5 and ~ 6 for 
IGR J16318-4848 and GX 301-2, respectively. These are very similar to the values obtained 
in previous studies for the total optical extinctions of the both sources. This suggests that a 
substantial portion, if not all, of the silicate features seen in Figures [1] and [2] are associated 
with the optical extinction of IGR J16318-4848 and GX 301-2. 



4. High-Resolution Spectrum: Line Emission 



Figures [3] and m present high-resolution spectrum of IGR J16318-4848 and GX 301-2, 
revealing numerous emission lines. (Note that the background emission is A^^OT subtracted.) 
The majority of them are H I lines, while the rest consists of metallic forbidden lines, pure 
rotational H2 lines, and 11.3 /im PAH emission feature. Table [T] lists the intensities of some 
of the bright H I lines identified in both sources with high (> 3) signal-to-noise ratio. The 
line intensities and the signal-to-noise ratio were estimated with Gaussian profile fits. The 
intensities of the H I lines are consis tent with gas at T ~ 10^ K and with A*"; 



10^ cm 



m 



the Menzel Case B state, following iFilliatre fc Chatvl (120041). altho ugh the intensity ratios 
are not very sensitive to those parameters (IHummer fc Storeyl 119871 ). 



Table [2] lists the parameters of forbidden ionic lines and H2 lines. For the forbidden lines, 
[Ne II], [Ne III], [S III], and [Si II] were detected in both sources, while [Ni II] and [Fe II] 
were only in IGR J16318-4848; [Fe III] is only in GX 301-2. The H2 lines were detected 
in both sources. The intensity ratio of the [S III] lines of 18.7 and 33.5 fim is diagnostic 
of the number density in the temperature range of 5000-20000 K. The extinction-corrected 
[S III] line surface brightness ratio of 18.7/ 33.5 urn is 1.21 for I GR J 16318-4848 and 1 68 fo r 
GX 301-2: according to the calculation of iHouck et al.l (119841 ) and [Alexander et al.l (119991 ). 
the ratios correspond to the electron number density iVe ~ 1000 cm~^ for both sources. 
On the other hand, the intensity ratio of the [Ne II] line (12.8 /im) and the [Ne III] line 
(15.6 /im) is sensitive to the hardness of t he ra diation, or temperature of the environment. 
If we use the calculation of iKunze et al.l (Il996l ). the extinction-corrected intensity ratio of 
[Ne II] /[Ne III] is consistent with ~ 4 x 10^ K for IGR J16318-4848 and ~ 4.5 x 10^ K for GX 
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301-2. The three H2 hnes detected in the both sources are pure rotational transition hues of 
H2 (0,0) S(0,l,2) which are usually excited by collisions. Assuming that the bottom two lines 
follow the distribution of the local thermodynamic equilibrium and are optically thin, the 
surface brightness ratios are consistent with the excitation temperature Tgx ~ 1000 K. The 
estimated H2 column densities are ~ 1 x 10^^ cm~^ (IGR J16318-4848) and 5 x 10^^ cm~^ 
(GX 301-2). This implies that the optical extinction associated with t he warm mo lecular 
gas is very small. Ay <C 1, based on the relation = 1.87 x 10^^ Ay ( lDrainell2003l ). 



5. Discussion and Conclusion 

Our SPITZER spectroscopic observations of IGR J16318-4848 and GX 301-2 have re- 
vealed for the first time the rich mid-IR environment of highly-obscured X-ray binaries. This 
includes two dust components with prominent silicate absorption, numerous H I recombina- 
tion lines, many forbidden ionic lines, and pure rotational H2 lines. Based on the observed 
spectra, we infer the following components for IGR J16318-4848 and GX 301-2: (1) hot (T 
> 700 K) and warm (T ~ 180 K) circumstellar dust; (2) ionized stellar winds responsible 
for the H I lines; (3) extended low-density ionized regions for the forbidden lines; and (4) 
photo-dissociated regions associated with the PAH, H2 and possibly the [Si II] line emission. 
For the forbidden lines, all the detected lines have relatively low ionization potentials like in 
starburst galaxies where the radiation is relatively soft (compared with active galactic nu- 
clei, for instance). This may indicate that the illumination of hard X-rays from the central 
compact X-ray source is not primarily responsible for the forbidden line emission. However, 
the inferred temperature for the radiation exciting the [Ne II] and [Ne III] lines is hotter 
than the stellar photospheres, so there can be some contribution from the compact object. 
Considering that [Ni II] and [Fe II] were detected only in IGR J16318-4848 while [Fe III] was 
only in GX 301-2, the radiation field of IGR J16318-4848 may be softer than GX 301-2, as 
demonstrated by the small temperature difference between the two sources (see §|4]). 

Perhaps the most natural explanation for the origin of the hot dust component relies on 
dust formation in the dense outflows from the early-type companions of IGR J16318-4848 and 
GX 301-2, as is the case in most sgB[e] stars (i.e., B-type supergiants with forbidden emission 
lines). However, the origin of the warm circumstellar dust component is very uncertain, and 
B[e] stars seldom show evidence for this type of warm dust. If both the hot and warm 
dust components have spherical shell geometres around the central star, then the associated 
optical extinctions are: Ay ~ 4 x 10"^ Qabs (Md/lO"*^ Mq) (Td/100 K)^ (Luv/10^^ ergs 
s~^), where and Td are the mass and temperature of the dust components, Qabs < 1 is 
the dust absorption coefficient, and Luv is the ultra-violet luminosity of the central star. 
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This gives the optical extinctions Ay <^ 1 for both the hot and warm dust components of 
IGR J16318-4848 and GX 301-2. (Here we use 1 x 10^^ ergs for Luv for both sources.) 
Therefore, under the assumption of spherical shell geometry, both the hot and warm dust 
components of IGR J16318-4848 and GX 301-2 contribute very little to the total optical 
extinction, suggesting that the hot and warm dust components are A/'OT strongly associated 
with the silicate absorption features. This also applies to the dust associated with the warm 
extended H2 gas since the optical extinction from this component is tiny (i.e., Av <^ 1; see 
§11]). What's then the origin of the silicate absorption features? If it is due to the foreground 
ISM, we would expect to see the CO2 ice feature around 15 /im, especially for IGR J16318- 
4848, based on th e intensity ratio between the sihcate absorption feature and ice feature 
found in the ISM (IKnez et al.ll2005l ). The absence of the ice feature in our spectra supports 
the interpretation that the silicate absorption features are probably not associated with the 
foreground ISM. One possibility may be the existence of an undisclosed colder (e.g., <C 100 K) 
circumstellar dust component which is responsible for the silicate absorption features and 
most of the optical extinction. We need further longer wavelength observatons to confirm this 
possibihty. Considering that the 9.7 /im silicate absorption feature represents the oxygen- 
rich material, the existence of the 9.7 micron silicate absorption feature may indicate that 
the origin of the potential colder dust component is related to the nucleosynthesis of the 
progenitors of IGR J16318-4848 and GX 301-2 (or their companions). 

The optical/near-IR companion of IGR J16318-4848 is a sgB[e] star. Such stars are 
known to have hot (T ~ 1000 K) circumstellar dust, and probably are evolving into Lumi- 
nous Blue Variables (LBVs) or Wolf-Rayet stars. Our mid-IR spectra of IGR J16318-4848 
and GX 301-2 are very si milar to that of the LBV P Cygni which shows many H I lines 
and forbidden ionic lines (ILamers et al.l Il996bl ) . The difference is that while the mid-IR 
continuum of P Cygni is due to the free-free emission in stellar winds, the main mid-IR 
emission of IGR J16318-4848 and GX 301-2 is from multiple dust continua. (The SEDs 
of both sources observed here are inconsistent with the free-free emission, as mentioned in 
Paper I.) However, we note that some LBV stars have also been observed to have thermal 
mid-IR dust emission ( ILamers et al.lll996al ). Based on the fact that B[e] stars seldom show 
mid-IR forbidden line emission, the B-type supergiant (or hypergiant) companions of IGR 
J16318-4848 and GX 301-2 may be in the evolutionary track to LBVs, implying that the 
extremely high obscuration seen in some massive X-ray binaries may be a phenomenon as- 
sociated with the evolutionary phase. This scenario is also consistent with the fact that 
the highly-obscured massive X-ray binary Cygnus X-3 has a Wolf-Rayet star companio n, 
together with the circumstellar dust emission of T ~ 250 K (IKoch-Miramond et al.ll2002l ). 



This work is based on observations made with the Spitzer Space Telescope, which is 
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Table 1. Observed Intensities of H I Transitions 



Transition (/xm) 


ICR 


GX 


Transition (/xm) 


IGR 


GX 




J16318 


301-2 




J16318 


301-2 


17-9(10.26) 


4.39 


2.44 


12-8(10.50) 


11.57 


11.40 


16-9(10.80) 


5.60 


2.20 


9-7(11.31) 


25.43 


13.50 


22-10(11.49) 


1.74 


0.78 


15-9(11.54) 


6.48 


1.87 


20-10(12.16) 


2.84 


1.35 


7-6(12.37)'^ 


72.02 


42.12 


18-10(13.19) 


5.15 


1.70 


13-9(14.18) 


10.61 


3.37 


23-11(14.30) 


2.00 


0.93 


22-11(14.71) 


5.02 


0.87 


16-10(14.96) 


5.41 


1.75 


20-11(15.82) 


2.79 


1.33 


10-8(16.21) 


14.66 


7.32 


15-10(16.41) 


3.83 


1.83 


19-11(16.59) 


2.15 


0.40 


12-9(16.88) 


7.35 


3.17 


14-10(18.62) 


4.04 


1.38 


16-11(20.92) 


2.69 


1.14 


9-8(27.80) 


18.71 


8.86 









Note. — The values are in 10 ergs s ^ cm ^ and are not corrected 
for extinction. 

^The H I 11-8 transitions may be blended. 
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Table 2. Observed Intensities of Other Lines 



Line(//m) 


ICR 


GX 


Line(//m) 


ICR 


GX 




J16318 


301-2 




J16318 


301-2 


[Ni II] (10.68) 


10.25 




[Ni II] (12.73) 


13.33 




[Ni II] (18.24) 


1.69 










[Ne II] (12.81) 


21.44 


23.41 


[Ne III] (15.56) 


2.60 


15.67 


[Fell] (17.94) 


22.54 




[Fell] (24.52) 


7.33 




[Fell] (25.99) 


70.93 




[Fell] (35.35) 


15.46 




[Fe III] (22.93) 




1.99 








[S III] (18.71) 


7.73 


8.02 


[S III] (33.48) 


40.14 


25.01 


[Si II] (34.82) 


177.20 


33.33 








H2(0,0)S(2)(12.28) 


2.88 


1.44 


H2(0,0)S(1)(17.03) 


4.66 


2.05 


H2(0,0)S(0)(28.22) 


10.08 


4.24 









Note. — The values are in 10 ergs s ^ cm ^ and are not corrected for 
extinction. 
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Fig. 1. — Main Panel: The background-subtracted IRS low-resolution spectrum of IGR 

J16318-4848. The solid line represents the observed spectrum, while the dotted line shows 
the best-fit spectrum combining the hot dust (dashed line), the stellar (dottcd-dashcd line), 
and the warm dust (dot-dashed line) components. The SED fit was conducted excluding the 
wavelengths associated with the silicate absorption features indicated by the arrows around 
9.7 and 18 /im. Bottom Panel: The residuals of the SED fits. Note that for the residual 
the fluxes of the wavelengths associated with the silicate absorption features were calculated 
with the best-fit parameters (see text). The unit of the residual is Jy. 
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Fig. 2. — Same as Figure [H but for GX 301-2. The silicate absorption features are not as 
prominent as in IGR J16318-4848, but are still visible, especially around ~9.7 fim. 
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Fig. 3. — The IRS high-resolution spectrum of IGR J16318-4848. Note that the background 
is A^Or subtracted. Some strong hues are marked. 
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Fig. 4.— Same as Figure El but for GX 301-2. 



